Unconditional spin squeezing via measurement-based quantum feedback 
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We demonstrate an unconditional spin squeezing in a collective spin system via feedback control 
based on quantum non-demolition measurement (QNDM) . We perform shot-noise limited collective 
spin measurements on an ensemble of 3.7 x 10 5 laser-cooled 171 Yb atoms in their spin-1/2 ground 
states. Correlation between two sequential QNDMs indicates generation of conditional spin squeezed 
states with — 0.801q'i2 dB of quantum noise reduction. Our feedback control successfully converts 
the conditional squeezing into the unconditional squeezing without significant loss of the squeezing 
level. 



Feedback is an essential building block of classical con- 
trol procedure. One can stabilize dynamical behaviour 
of the target system by sensing its state and manipulat- 
ing the system depending on the sensed-outcomes. Re- 
cently, the real-time stabilization of photonic quantum 
system [1; has been also demonstrated by using quantum 
non-demolition measurements (QNDMs) and the feed- 
back control. Such an active control of quantum system 
is a significant step towards realization of advanced quan- 
tum information processing [2H1]- 

In addition to photonic systems, atomic quantum sys- 
tems are also very attractive from a viewpoint of quan- 
tum feedback. The measurement-based quantum feed- 
back of an atomic spin ensemble, in particular, enables 
us to realize important applications such as the quantum 
memory for continuous- variable system [5H7] and the en- 
hancement of the squeezing of the squeezed spin state [5] 
which would improve high precision spectroscopy such as 
an optical lattice clock [9] or an atomic magneto-meter 
[T01 - H2] . While there are several experiments that achieve 
spin-squeezing of an atomic ensemble by using QNDM 
[13H18j . there has been no report on the measurement- 
based quantum feedback control of a squeezed spin state. 

Here we have achieved unconditional squeezing of the 
atomic quantum fluctuation via successful feedback con- 
trol of an atomic collective spin state by using QNDMs 
and coherent manipulations |19) . Our feedback control 
successfully converts the conditional squeezing into the 
unconditional squeezing without significant loss of the 
squeezing level. Our result will open the door towards 
emerging field of quantum information dynamics with 
feedback 

We note that our measurement-based feedback is not 
'coherent quantum feedback' [23H2Z], with which one 
controls the dynamics of the target system without mea- 
surement device. Particularly interesting case is the 
cavity-feedback scheme [26l |27], where an impressive 
squeezing level for atomic spins is achieved. Compared 



with the coherent feedback, our approach does not re- 
quire any purpose-built interaction. The power of the 
quantum measurement enables us to design various (es- 
sentially unlimited variation of) time evolution of quan- 
tum system. 

A target quantum system in this work is an atomic 
collective spin described by a collective spin vector J = 
(Jx, Jy, Jz)- The uncertainties of spin components, SJ y 
and SJ Z , are jointly limited by the uncertainty relation 
SJ y ■ SJ Z > I (J x ) |/2 due to the commutation relation 
[Jy, Jz] — iJx- In our case, the initial state is given by an 
iYi-partite product state \<p) = \*f x ) (g) \-f x ) <g) • • • <g) |t a ), 
where |ts) is an eigenstate of the x component of a single- 
particle operator of spin-1/2. The state 1^) is an eigen- 
state of J x with the eigenvalue of J = Y4/2, and its 
standard deviations of the transverse spin components 
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FIG. 1. (Color online) Schematic illustration of the controlled 
dynamics. The atomic spin fluctuation is shown as a circle 
on the surface of the collective Bloch spheres. The shadows 
projected on the screens represent the total spin fluctuations. 
The quantum non-demolition measurement (QNDM) brings 
the initial atomic state into the ensemble of various squeezed 
spin states, and the quantum feedback (QFB) displaces the 
squeezed spin states so as to suppress the total spin fluctua- 
tion along the z-axis. 
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are obtained as 5J y = SJ Z = J/2. It therefore shows 
the minimum uncertainty SJ y ■ SJ Z = J/2 equally dis- 
tributed over any two orthogonal components perpen- 
dicular to the x direction, and \<p) is referred to as the 
coherent spin state. 

A QNDM on the atomic collective spin by using a po- 
larized light field as a probe system [551 HH] brings the 
coherent spin state into other quantum states with re- 
duced quantum uncertainty [16 , which is schematically 
depicted in Fig. [I] The polarization degrees of freedom 
of light field is described by the photonic Stokes vector S 
whose components are defined as S x — (aia y — a x a x )/2, 
S y = (a y a x + a x a y )/2 and S z = (a y a x — a x a y )/2i, where 
a x and a y are annihilation operators of light field linearly 
polarized along the corresponding directions. The quan- 
tum non-demolition interaction between these two sys- 
tems is implemented by the Faraday rotation, by which 
the axis of the polarization of light is rotated by an an- 
gle being proportional to the measured atomic collective 
spin. The interaction is described by a unitary operator 
Uf = cxp (— ixf J Z S Z ), where \f represents the Faraday 
rotating angle per unit spin angular momentum. The 
QNDM induces a type of inter-particle correlation or en- 
tanglement in the collective spin system, and reduces the 
fluctuation of J z below J J/2. Such a state is known 
as a squeezed spin state [3U]. However, the obtained 
squeezed spin states is randomly distributed depending 
on the measurement outcomes over the range of the ini- 
tial quantum fluctuation of \<p) as in Fig. [TJ the total 
fluctuation along the z-axis just after the QNDM is the 
same as that of the initial coherent state. The goal of our 
work in this paper is to suppress this random fluctuation 
via quantum feedback (QFB) OH]. 

For the experimental implementation of the QFB, we 
use an ensemble of ytterbium ( 171 Yb) atoms as the tar- 
get system. The 171 Yb atoms have only a nuclear spin- 
1/2 in their ground state. The system therefore consti- 
tutes the simplest 2-level energy structure and is robust 
against decoherence due to residual magnetic field owing 
to the small nuclear magnetic moment. There are no ef- 
fects of an additional rank-2 tensor term which appears 
in the case of atoms with larger spin, such as usually 
employed alkaline atoms. Figure [2] shows the schematic 
of experimental setup. We firstly create the coherent 
spin state by a 10-//S long optical pumping pulse of cir- 
cularly polarized light propagating along the x axis and 
tuned to the 1 S Q <-> 1 Pi(F' = 1/2) transition. Then we 
measure the rotation angle of linearly polarized probe 
pulses S y which get the information of the atomic col- 
lective state via the Faraday rotating interaction. The 
probe pulses are tuned to the 1 So «-> 1 Pi(F / = 1/2) 
transition with A = — 2tt x 160 MHz detuning, and the 
mean photon number per the pulse is Nl = 1.3(1) x 10 6 . 
The probe is horizontally polarized before the interac- 
tion (S = | (S x ) | = N L /2). The Faraday rotating angle 
with the atoms spin-polarized along the probe direction is 




FIG. 2. (Color online) Apparatus to implement the quantum 
feedback control. AOM, Acoustic- Optical modulator; PBS, 
Polarization beam splitter; PDs, Photo detectors. Horizon- 
tally polarized probe pulses illuminate the ensemble of cold 
171 Yb atoms which are spin-polarized along the x direction by 
an optical pumping beam. The polarization and Sy of 
the pulses are measured by a polarimeter which consists of a 
PBS and two PDs. A time interval between these two pulses 
is 12.5 fis. The measurement outcomes are fed into the logger 
and feedback controller which determines the pulse duration 
T(Sy ' 1 ) of the circularly polarized control pulse. 



0.16(2) rad, which corresponds to an effective atom num- 
ber of Na = 3.7(4) x 10 5 . 

The QFB in this work proceeds in the following way; 
the output of the polarization detector S y °^ is fed into 
both of data logger and feedback controller, then the out- 
put is looped back into the atomic system. The feedback 
operation is implemented by a fictitious magnetic field 
[31] which is generated by a circularly polarized 'control 
pulse' propagating along the y-axis. After the QFB, we 
measure the spin state by applying the second QNDM 
pulse Sy 1 ^ . A time interval between the first and the sec- 
ond QNDM is 12.5 fis which is sufficiently shorter than 
the typical decoherence time (several hundred microsec- 
onds due to ballistic expansion of cold atomic cloud) . 

In the conventional normalized quadratures for light 
and atoms polarized along the a;-axis (Xl,Pl) = 
(S v ,S z )/y/S and (Xa,Pa) = {Jy,Jz)/\ /f J, we can work 
with the standard commutation relation [Xl,Pl] — 
[Xa,Pa] = i and uncertainty relations V(Xl)V(Pl) > 
1/4 and V(Xa)V(Pa) > 1/4 from the commutation re- 
lation, where V(-) means the variance. By using these 
operators, the Faraday rotating unitary operator Uf can 
be rewritten as Uf — e~ lX}JzSz = g ~ ikP aPl w ith a di- 
mensionless parameter k = x/V SJ, which is calculated 
as PU 

2CT r v / 5J A _ s 

K ~ 37rwg(l + s) A 2 + (r/2) 2 ' S ~ l + (2A/r) 2 ' 

where T — 2n x 29 MHz is the natural linewidth, <r — 
7.6 x 10 _14 m 2 is the photon-scattering cross section in 
171 Yb atom, the beam waist wq = 40 /im, and the sat- 
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FIG. 3. (Color online) Normalized variances of the sum (<$+, 

A) and the difference (5_, v) of two outcomes S v and 
as a function of feedback gain (the error bars are the statis- 
tical standard deviations Ha). The solid curves are the cor- 
responding theoretical predictions from the calculation scaled 
by a single proportional factor to link the measured intensity 
of control light pulse to feedback gain. The insets show mea- 
sured joint frequency distributions of two sequential QNDM 
outcomes (horizontal axis) and Sy 1 ^ (vertical axis) with 
solid curves indicating 3a radii. 

uration parameter sq — 7.2. In our case, k is estimated 
as 0.59. For these parameters, the damping coefficient 
[32] of light {e£) and atoms (e^) can be estimated as 
(e L ,e A ) = (0.042,0.15). 

The control pulse, which is tuned to the 1 So *H> 
3 P 1 (F' = 1/2) transition with -2tt x 70 MHz detuning, 
induces the phase shift and rotates the collective atomic 
spin vector in the x-z plane. The rotation angle is pro- 
portional to both of the time duration and the intensity 
of the control pulse. Our feedback controller, which is 
implemented by a field-programmable gate array, deter- 
mines the time duration depending on the shot-by-shot 
measurement outcome. The intensity of the control pulse 
is an additional parameter of our feedback, and is used 
as 'feedback gain'. The negative sign of the gain assigned 
for the correct feedback direction to successfully reduce 
the fluctuation. The direction of the feedback can be 
changed by switching the sense of the circular polariza- 
tion of the control pulse. 

Figure [3] shows the observed correlations between the 
two sequential QNDM outcomes modified by the feed- 
back control. We show the normalized variances of 
the sum and the difference of two outcomes as a func- 
tion of the feedback gain, which are defined as <5± = 

V{xf ] ± Xi X) )/2 = V{S ( y ] ± S y 1) )/{2S). Each point is 
calculated based on 10, 000 independent sets of the out- 
comes with the corresponding feedback gain. The mea- 
sured joint frequency distributions of two outcomes S y °^ 
(horizontal axis) and Sy (vertical axis) are also plotted 



at some feedback gains. The normalized variances S± 
correspond to the width of the distribution along diago- 
nal lines of the individual frame of the insets. The value 
of 0.5 of 6± corresponds to the light shot noise. Here we 
also show the actual data which is obtained without the 
atoms, and we can see that the normalized variances is 
close to the expected value. The excess spin noises are 
clearly observed for the data with the atoms, regardless of 
the feedback gain. The data with zero feedback gain in- 
dicates generation of the conditional squeezed spin states 
[HO |3I] ■ We have confirmed that the observed noise is 
of quantum origin by several ways |33j . A comparison 
between 5i and St highlights the effect of the feedback 
control. As schematically shown in Fig. [T] the atomic 
spin states are expected to be a certain squeezed spin 
state via the optimal QFB. It means that the atomic 
state after the optimal QFB does not depend on the first 
outcome, the correlation 5\ — o 2 . is therefore expected to 
be close to zero at the optimal feedback gain g = g op t- In 
other words, our QFB reduces the quantum fluctuations 
at the expense of the correlation or the information. If 
\g\ < | <? pt | which means that the rotation angle of the 
mean spin is smaller than the optimal one, the correla- 
tion does not reach zero but is still positive. On the other 
hand, if \g\ > |g pt|, the correlation becomes negative one 
due to the overshooting of the rotation. Figure [3] shows 
a transition of the correlation from positive to negative 
with increasing the intensity of control pulse |<?| at g < 0. 
This is a direct evidence that the atomic spin states can 
be manipulated depending on the shot-by-shot measure- 
ment outcomes. The results are in good agreement with 
our theoretical model [M] (solid lines in Fig. [3]). 

In order to know the available maximum squeezing 
level of the unconditional spin squeezing in our case, we 
evaluate a conditional squeezing parameter. From the 
positive correlation between two QNDM outcomes in the 
case of no feedback, we define 

2 _ ^cond — V{S y o) 

where V con d = min V(Sy + g c Sl ' ) ), and V^(sIq^) cor- 
gc y 

responds to the variance of the second QNDM in the 
case of no atoms. The conditional variance V^ond means 
the available minimum variance from the information ob- 
tained by the first probe. Therefore, C^ond represents how 
much the atomic spin fluctuations can be suppressed by 
the correlation between Sy ^ and Sy with virtual and 
ideal QFB. The optimal value of g c for the minimiza- 
tion is given by ~cov(S y °\ Sy 1 ^)/V(S y ^), where the cov 
means the covariance. The conditional squeezing param- 
eter is estimated as ^ ond = 0.83 ± 0.02 = -0.801^ dB, 
which indicates the successful conditional spin squeezing. 
The experimental setup for the conditional squeezing is 
basically the same as that of our previous works [35 . 
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FIG. 4. (Color online) (a) Unconditional squeezing parame- 
ters £u Ilc as a function of feedback gain, (b) For comparison, 
the first QNDM outcomes S y are used in the calculation 
of the numerator of the £u nc , which corresponds to that of 
no feedback (£; of ). The error bars represent the statistical 
error ±1<t. The solid curves are the theoretical predictions. 
The filled region shows the conditional squeezing parameter 

dud = — 0.80 



+ 0.11 
0.12 



dB to indicate the optimal squeezing. 



With the QFB, a measure of the unconditional spin- 
squeezing is simply the variance of the second outcome 
Sy' . Here we define the unconditional squeezing param- 
eter £^ nc which characterizes the effect of the QFB as 

c2 V(SP)-V(S$) 



y(4 ref) ) 



The important feature of the QFB is that we can ma- 
nipulate the dynamics of the target quantum system via 
consecutive feedback. While the achieved squeezing level 
and the successful feedback steps in this study are lim- 
ited because of the absorption of photons in probe pulses 
and the ballistic expansion of the atomic cloud, these 
problems are not of the fundamental origin and can be 
avoided by realistic improvement such as using an opti- 
cal trap system to increase the optical density of atomic 
cloud. We believe that our results pave the way for more 
advanced QFB [3J 21 [TU] of the collective spin system. 

In summary, we have experimentally demonstrated the 
unconditional spin-squeezing via quantum feedback con- 
trol which converts the initial coherent spin state to the 
squeezed spin state in an unconditional manner. An un- 
conditional squeezed spin state can be directly used as 
a resource for quantum-enhanced metrology. We note 
that our feedback controller can be regarded as a quan- 
tum analog of Maxwell's demon [2T)lJ2"2l 155] , While the 
original classical demon only suppresses thermal fluctu- 
ations, the quantum demon can also suppress quantum 
fluctuations. Moreover, the inter-atomic entanglement 
in squeezed spin states in our experiment can be a useful 
resource for quantum heat engines [35] . 
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work is supported by the Grant-in- Aid for Scientific Re- 
search of JSPS (No. 18204035, 21102005C01 (Quantum 
Cybernetics)), GCOE Program "The Next Generation 
of Physics, Spun from Universality and Emergence" from 
MEXT of Japan, and World-Leading Innovative R&D on 
Science and Technology (FIRST). 



where Sy ie ® and S y T Q^ are the outcomes of the polariza- 
tion measurements with and without atoms in the case 
of no feedback. They are measured just before the main 
sequence 



j , although it is still possible to use S y ^ and 



S y °Q instead of S y luI> and SyQ 1 ' , respectively. The uncon- 
ditional squeezing parameter £„ nc represents how much 
the atomic spin fluctuations are suppressed via the actual 
QFB. 

Figure Qa) shows the unconditional squeezing param- 
eter as a function of the feedback gain. While no un- 
conditional squeezing is observed for the case of no QFB 
(b), we clearly observe the maximum reduction of the 
spin fluctuation in the vicinity of the feedback gain where 
the correlation 6+ — is close to zero (see also Fig. |3j) , 
as expected. This demonstrates the unconditional spin- 
squeezing as a result of the successful QFB. In addi- 
tion, we have successfully observed the features such as 
maximum (and near-optimal) squeezing at some negative 
gain, no squeezing at zero and too much negative gain, 
and also the increase of the noise at positive gain. The 
overall features of the observation in Fig. [4] are in good 
agreement with the theoretical model [37]. 
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EXPERIMENTAL SETUP 

In this study, we use a cold 171 Yb atomic cloud as 

the target system of our quantum feedback control. The 

whole setup and timing sequence are shown in Fig. [5] 

The measurements are performed within the "MOT-off" 

periods, where the magnetic and light fields responsible 

for cooling and trapping are switched off. During the 

period, we firstly create the coherent spin state by a 10-jus 

long optical pumping pulse of circularly polarized light 

propagating along the x axis and tuned to the 'So f> 

1 Pi(F' = 1/2) transition. Then we measure the rotation 

angle of linearly polarized probe pulses which get the 

information of the atomic collective state via the Faraday 

rotating interaction. The probe pulses are tuned to the 

% O 1 Pi(F' = 1/2) transition with -2ir x 160MHz 

detuning as shown in Fig. [5jb) . 

Figure pic) shows the main-sequence of the experi- 
(•) 

ment. Here S y denotes the outcome which is obtained 
within the main-sequence. Within the 'Measurement 
& Feedback' region depicted by the dashed square in 
Fig. [5|c) , we perform quantum feedback experiment by 
utilizing the control pulse, which is circularly polarized 
and tuned to the 1 So <-> 3 Pi(F' = 1/2) transition with 
— 2tt x 70 MHz detuning. The time- width T of the control 
pulse is determined by the foregoing outcome Sy \ and 
the function is implemented by the field-programmable 
gate array as the 'Feedback Controller' and by the 'AOM' 
(Acoustic-Optical Modulator) as schematically shown in 
Fig. [5ja). In this system, we can control the feed- 
back gain g by changing the intensity and the polariza- 
tion of the control pulse. The main sequence was re- 
peated 10, 000 times with a certain value of the feedback 
gain. This is a unit experimental run of the experiment, 
and is also repeated with various values of the feedback 
gain. After the measurement and control period, the new 
atomic cloud is formed for the next sequence at a rate of 
about 100 Hz. 

We also measure the light shot noise by using the same 
setup but with no atoms just before every experimental 
run; this is the sub-sequence of the experiment and is 
shown in Fig. [5jc'). In the sub-sequence, we turn off 
the magnetic field not only within the 'MOT-off' period, 
but also within the 9.7-ms long loading period so as not 
to load the atoms. The outcomes obtained in the sub- 
sequence, which are denoted as S y Q , are used to estimate 
the amount of conditional/unconditional squeezing. 



Performance of QNDM 

Figure|6ji)-(iii) show the correlation plots of the two se- 
quential quantum non-demolition measurement (QNDM) 
outcomes with (i) no atoms, (ii) independent coherent 
spin states, and (iii) identical collective spin states in 
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FIG. 5. Schematic illustration of the experimental setup, (a) Experimental apparatus, (b) Energy level structure of 
171 Yb and associated laser frequencies, (c) Main-sequence of the experiment, (c') Sub-sequence of the experiment. The details 
are described in the text. 





(i) without atoms 


(ii) independent spin states 


(iii) identical spin states 




3.21(1) x 10 5 


4.43(2) x 10 5 


4.44(2) x 10 s 


V(S^) 


3.22(1) x 10 5 


4.44(2) x 10 5 


4.43(2) x 10 s 


A% 


3.29(1) x 10 5 


4.73(2) x 10 5 


5.39(2) x 10 s 


Al 


3.15(1) x 10 5 


4.15(1) x 10 5 


3.49(1) x 10 5 


Scond 




0.98(2) 


0.83(2) 



TABLE I. Variance and conditional squeezing corresponding to Figj6] (i)-(iii). All values are calculated from 150,000 
polarimeter outcomes, which are converted to the corresponding photon number differences. The outcomes are collected 
simultaneously with the quantum feedback experiments. The numbers in parentheses denote la statistical errors. 
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the case of no feedback . In our experiments, the pairs 
of outcomes Sy ^ and Sy 1 ^ corresponding to the case (ii) 
and (iii) are obtained simultaneously with that for quan- 
tum feedback experiments as shown in Fig. [5jc); the 
pair of Sy ^ and Sy corresponds to the case (ii), and 
the pair of Sy and Sy corresponds to the case (iii). 
Here we confirm that the observed variances V(Si 0) ) and 
V^Sy 1 ^) are consistent with the quantum noise limited 
values. Each probe pulse of the QNDMs has on av- 
erage Nl — 1.3(1) x 10 6 photons. From this value, 
we can estimate the expected variance in the case of 
no atoms (i) as V(S y ) = N L /4 = 3.3(3) x 10 5 , which 
is consistent with the observations in Table I (i). The 
atom number in the probe region is determined from the 
measurement of the Faraday rotation angle 6 with the 
sample spin-polarized along the probe direction. The 
Faraday rotating angle per unit spin angular momentum 
Xf links the effective number of atoms with the Fara- 
day rotating angle as 9 — A typical value of 
6 = 0.16(2) rad for the rotation angle corresponds to the 
effective atom number of 3.7(4) x 10 5 . We can estimate 
the expected quantum-noise-limited variance of the po- 
larimeter outcomes with the sample spin-polarized per- 
pendicular to the propagation axis of the probe pulse as 
V{S y ) = N L /4+( X fN L /2) 2 -N A /4 = 4.4(4) x 10 5 , which 
is consistent with the experimentally obtained variances 
V(S y 0) ) and V(S y 1} ) = 4.4 x 10 5 listed in Table I(ii) and 
(iii). 

Here the correlation between two QNDM outcomes 
Sy ^ and Sy 1 ^ is clearly seen in Fig.^iii). Quantum noise 
reduction is now only possible in a conditional way via 
the correlation. The standard deviation A± which is de- 
fined as Aj_ = V(Sy ^ ± Sy 1 * > )/2 characterizes the corre- 
lation, and corresponds to the width of the distribution 
along diagonal lines of the individual frame of Fig. [6j The 
obtained Aj_ is summarized in Table I with statistical er- 
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-2-1012 -2-1012 -2-1012 

(xlO 3 ) (xlO 3 ) (xlO 3 ) 

Sf> [photons] Sf> [photons] S™ [photons] 



FIG. 6. Joint distributions of the two sequential 
QNDM-outcomes S y ^ and Sy without any feedback 

control. The measurement is performed with (i) no atoms, 
(ii) independent coherent spin states prepared by an addi- 
tional pumping pulse between two probe pulses, and (iii) iden- 
tical collective spin states. Solid curves indicate 3a radii and 
dashed curves in (ii) and (iii) are the same as the solid curve 
in (i). 



rors. One can clearly see a reduction of A_ as well as 
an enhancement of A + in the case of identical spin state 
(iii), compared with A± of independent spin state (ii) 

and also the respective V(S y 0S> ) and V(S y lS> ). 

In order to check that the observed noise is of quantum 
origin, we also evaluate the conditional squeezing param- 
eter with applying an additional pumping pulse during 
the two probe pulses Sy ^ and Sy . The pumping pulse 
resets the collective atomic spin state and thus the corre- 
lation due to the quantum fluctuation is expected to dis- 
appear as a result of the re-initialization. Consequently, 
the correlation, if observed even with the application of 
the pumping pulse, should indicate the contribution from 
the classical noise. In our experiment, we minimize the 
possible residual correlation by adjusting an incident an- 
gle of the pumping pulse, which results in the condi- 
tional "squeezing" parameter of ^ond = ± 0.02 = 
—0.068 ±0.10 dB for the independent coherent spin state 
shown in Fig. [6]Jii) . This result confirms that the effect 
of the classical noise is too small to mask the quantum 
fluctuations and our QNDM has sufficient sensitivity to 
access the quantum fluctuations. Possible origin of the 
classical noise is the residual magnetic field and the probe 
pointing fluctuation. The slight reduction of the condi- 
tional squeezing level compared with our previous work 
|16) comes from the slight decrease of the atom num- 
ber, and has nothing to do with the performance of the 
measurement-based quantum feedback control. We note 
that an increase of the conjugate quadrature noise due 
to the QNDM was observed in our previous work using 
essentially the same experimental setup [31 j . 

Theoretical model 
QND interaction with scattering loss 

After the probe pulse passes through the atomic en- 
semble, the QND interaction under the influence of the 
scattering loss changes the canonical operators Xl and 
Pa [22]: 

X' L = VT^l{X l + kP a ) + ^£E L , (1) 
P'a = VT^P A + V^E A , 

where the symbols with (without) a prime denote the 
operators after (before) the interaction, El and Ea are 
the standard vacuum operators with variance 1/2, and cl 
and €a are the damping coefficients. The vacuum oper- 
ators represent the contributions from unpolarized com- 
ponents due to the scattering. There are two essential 
assumptions to derive Eq. ([!]): (i) one-dimensional model 
is appropriate, in other words, the interaction volume is 
of a pencil shape with Fresnel number F ~ 1 and (ii) 
the excited atoms decay into either of the two spin states 
in the ground state. The scattered photons are emitted 



into the same spatial mode with the probe pulse because 
of the former assumption, and thus are fed into the po- 
larimeter. This effect is represented by the contribution 
of El . Similarly, all atoms which absorb a photon in the 
probe pulse come back to either of the two spin states in 
the ground state because of the latter assumption, there- 
fore the atoms are not lost but partially depolarized due 
to the photon absorption-and-scattering process. The ef- 
fect is represented by the contribution of Ea- 

In our experiment, the Fresnel number F ~ 10 is much 
larger than unity and there are no other ground states; 
i.e., the latter (ii) is appropriate, but the former (i) is 
not. The QND interaction in our experiment is therefore 
described by the following transformation: 

X' L = V\^l{X l + kP a ), (2) 
P'a = VT^Pa + V^E A . 

Uncorrelated noise 

Variance of the sum (5% = V{X L a) + X L 1] )/2) and 

the difference (51 = V(X { L 0) - X^)/2) characterize the 
correlation between the two sequential outcomes. For the 
ideal case (e A = =0) without feedback (g = 0), the 
difference X L °^ — becomes X L q — X^q and this is 
equal to the case without atoms (k = 0). The variance 
51 without feedback (g — 0) can also be calculated by 
using the Eq.([2]) as 

5 2 _(k, e A ,e L ) = (1 + k 2 {1 - VT^a)) , 

and 51 (n,e A ,e L ) < V(X$,) = 1/2 holds over a 
wide range of the parameters around (k,ca,^l) = 
(0.59,0.15,0.042), e.g., 51 (0.59,0.15,0.042) = 0.49. 
However, as shown in Table I (i) and (iii), the experi- 
mentally observed Al = 3.49(1) x 10 5 is slightly larger 
than the variance of the light shot noise 3.21(1) x 10 5 
which is measured without atoms. The results sug- 
gest that the measurement outcomes include some ad- 
ditional uncorrelated noise, which cannot be cancelled 
out by taking the difference between two sequential out- 
comes similarly to the light shot noise. The amount 
of the additional noise can be estimated as V(El) = 
(3.49(1) x 10 5 - 3.21(1) x 10 5 )/2 = 0.14(1) x 10 5 . We 
emphasize that the kind of noise cannot be reduced by 
using the feedback in any way, and it just decreases the 
amount of the noise reduction. We can effectively take 
into account the presence of uncorrelated noise by adding 
the uncorrelated noise terms in Eq.Q as follows: 

X' L = JT^V L (X L + kP a ) + J7 l E l , (3) 
P' A = VI^Pa + V^aEa. 

where a pair of e' L and El (with variance 1/2) repre- 
sents the contribution form the additional uncorrelated 



noise. The term yJe^E^ includes not only the scatter- 
ing photons as the standard vacuum operator in Eq.Q 
or Ref. 32J, but also the contributions from all uncorre- 
lated noise in the presence of atoms. We can estimate 
e' L = 0.098 from the modified model Eq.([3]) with the ex- 
perimental results Al = 3.49(1) x 10 5 and the variance 
of the light shot noise 3.21(1) x 10 5 . 

Quantum feedback 

Adding the feedback term to Eq. ^ yields the following 
formula: 

X' L = VT^T L {X L + kP a ) + ^ L E L , (4) 

P' A = VT^Pa + yfT A E A 

+ g-(l-e A )(X' L -x ), 

where the third term of P' A represents the effect of feed- 
back. The pre-arranged target value xq has set to be 
about 3 • yJV(X' L ) in the experiment. The absorption 
makes effective feedback gain to be 1 — ea times smaller 
than that in the case of no scattering loss. 

We have also done the multi-cycle feedback in the ex- 
periment, as shown in the timing sequence (Fig. [5jc)). 
We can express the multi-cycle feedback by using the 
following recurrence formula which corresponds to the 
generalized form of Eq. ^ as 

X? = y/T=I£{X® + K Pf) + y/7 L E®, (5) 

+ g-(l-eA)(xt 1) -x )-F[xt 1 \x ], 

r n f 1 a < b 
where F\a,b\ = < 

11 [0 a > b 

Here the variances are normalized as V(X^) — 
V(P A 0) ) = V{E { 1 ] ) = V{Ef) = 1/2, and the work- 
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-300 -200 -100 100 
Feedback gain (a.u.) 

FIG. 7. Two-cycle feedback results as a function of 
feedback gain. The error bars represent the statistical er- 
ror ±lcr and the solid curve is the theoretical predictions. 
The filled region shows the conditional squeezing parameter 
£cond ~ —0.80^0 '12 dB to indicate the optimal squeezing with 
1-cycle feedback (Fig. 4(a)). 
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ing range of our feedback is limited as represented by 
the function F[Xj? 1 \xo]. We define the conditional 
squeezing parameter with the i-cycle feedback as £f = 

(K(4Vi/2)/(n4 0) ) : i/2). 

Figure [7] shows the experimental result of £f , and it is 
again in good agreement with the theoretical prediction 
calculated by Eq.([5|. Note that the multi-cycle feedback 
is performed with the pre-fixed feedback gain, which is 
known as not the optimal one [19] . Multi-cycle feedback 
is expected to improve the squeezing and broaden the 
gain dependence. 



Feedback control with partially polarized spins 

In the theoretical model described above, we assume 
that the unpolarized component of atoms (t/caEa) is not 
changed by the feedback control. Here we consider the 
reliability of the assumption. First of all, it is likely that 
we cannot get any information of the unpolarized spins 
from the first outcome xf\ and thus the feedback con- 
trol, which depends on , only applies some random 
spin-rotation and does not contribute the noise reduction 
of the unpolarized spins. In addition, even if we could get 
the information of the unpolarized spins from X^ , the 
amount of the information is much smaller than the light 
shot noise because the upper bound of the signal to noise 
(signal from the unpolarized spins to the light shot noise) 
ratio can be estimated as k 2 ■ ea — 0.05. 

If we could obtain the information of the unpolarized 
spins with sufficient signal to noise ratio, the noise of the 
unpolarized spins can be reduced by the feedback. How- 
ever, the gain dependence with the unpolarized spins is 
supposed to be drastically different from that with po- 
larized spins. For the unpolarized spins, the positive 
and negative feedback equally reduce the spin fluctua- 
tion whereas the negativity of the feedback gain is es- 
sential for the polarized spins. This is because there 
are no special directions for the unpolarized spins ex- 
cept for the measurement axis, whereas the mean spin 
direction is well-defined for the polarized spins. In ad- 
dition, the amount of the noise reduction for the unpo- 
larized spins is negligibly small comparing with that for 
the polarized spins in the parameter range of the ex- 
periment. This can be qualitatively understood in the 
following way. In order to change the distribution of the 
z-component of angular momentum m z , we displace m z 
by rotating the spins in our setup. The amount of the 
displacement with small rotation angle j3 can be esti- 
mated as 5m z = j/3, where j is the total angular momen- 
tum. The polarized spin-1/2 system with N atoms can 
be represented by the number of (N + 1) state-vectors 
with j = N/2 = J while the unpolarized mixed state re- 
quires 2 N state- vectors whose total angular momentum is 
mostly smaller than J. For the noise reduction with the 



unpolarized spins, the rotation angle /3 therefore has to 
be larger than that with the polarized spins to achieve the 
same amount of the displacement. The mean total angu- 
lar momentum of the unpolarized mixed state can be es- 
timated as ~ 0(J{J%) + (Jf ) + ( Jf )) ~ o(VJ). 
In our feedback control, the rotation angle j3 is given by 
the measurement outcome m z and the feedback gain g as 
f3 = gm z . For the polarized spins, the feedback with the 
optimal gain g oy which gives the maximum noise sup- 
pression, is determined by /3 = g m z and J/3 — —m z , 
then the optimal feedback gain is roughly \g a \ ~ 0(1/ J). 
This is consistent with the experimental parameter in our 
feedback control. On the other hand the optimal gain 
for the unpolarized spin state, say g' Q , is much larger be- 
cause the larger rotation angle is required with the same 
outcome m z . Similarly, the optimal value is determined 
by the inverse of the mean total angular momentum as 
\g' \ ~ 0(l/y/~J). Since J ~ 10 5 in the experiment, the 
value of \g' a \ estimated from this argument is more than 
100 times larger than the value of \g a \. 

The asymmetric behaviour of polarized spins is clearly 
observed in Fig. 3 and 4. It confirms that the observed 
spin noise reduction arises from the feedback control of 
the polarized component, and that the most of the unpo- 
larized component is not changed within the parameter 
range of the feedback gain, which is around the optimal 
one for polarized spins \g \ <C \g' \. 

In the following part of this section, we analyze the ef- 
fect of the feedback control on the unpolarized spins more 
quantitatively. In order to see the effect of noise suppres- 
sion on the unpolarized spins, let us define the complete 
basis labelled by j, m z , and an additional quantum num- 
ber d. Using the basis m z ,d)}, we can express the 
unpolarized iV-spin state as 

Pn = I 2 n/2 n (6) 

N/2 j P Nj 
j=jo m x =—j d=l 

where jo = with the even N, and jo = 1/2 with the odd 
N. The magnitude j is either an integer or half-integer 
and is less than or equal to J — N/2. The quantum 
number d = 1, • • • , Pj^j distinguishes between the P/vj- 
degenerated subspaces |4"0] : 

N\(2j + 1) 

Nj (N/2-j)\(N/2+j + iy: {) 

The mixed state Pat is already given by the ensemble 
of the quantum states which are characterized by m z . 
Since m z corresponds to the first measurement outcomes 
4 0) or xf , we use the character m instead of m z for 
the sake of convenience. Note that we ignore the read- 
out noise here, but the approximation just emphasizes 
the gain-dependence for the unpolarized spins. We can 
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FIG. 8. Numerical calculation of the measurement-based feedback for completely mixed states, mo and mi 

correspond to Sy ^ and Sy 1 ' in the experiment (in the case that the read-out noise can be negligible), respectively, (a) Variances 
of the sum (6+) and the difference (S'L) calculated with pioo as a function of feedback gain, and corresponding joint probability 
distributions (insets), (b) Unconditional noise reduction after the feedback for pioo as a function of the feedback gain, which 
is essentially the same as the plot in Fig. 4(a). 



express the state after the feedback as 



function [41] D%> 0>mi (fi) = (j, mi, d|e jJ ^|j, m , d) as 



p(mo, mi) = Tv[Ml x M miP '] = Tr[M miP '} (8) 



N/2 



2 N 



3=]a 



! N/2 j P Nj 

;]v X! U (9 m o) \j, m o,d)(j,m ,d\ U\gm ), 



2 N 



j=jo m =-j d=l 



where g is the feedback gain, and U (/?) = exp (— i J y /3) 
represents the rotation about the y-axis with the ro- 
tation angle /3. Our interest is in the probability of 
getting the outcome mi, which corresponds to Sy 1 ^ 
or Xjp, by measuring the observable J z of the state 
p' . Corresponding measurement operator is given by 



M„ 



_ ^N/2 



Ej=j T,d=i \ j, ™>i,d){j, mi,d\. The joint prob- 
ability distribution of mi (after the feedback) and mo 
(before the feedback) can be calculated with Wigner D 



Probability distributions p(m ) and p(mx) can also be 
calculated as p(mi) — J2 ma P(' m 0i rn i) an d p( m o) — 
E mi P(m , mi ). 

We can estimate various variances V(-) by using the 
probability distributions. Figure |8ja) shows the numer- 
ically estimated variances of the sum (6+ = V(m + 
mi)/2J) and the difference (d^_ = V(m — m\)/2J) as a 
function of the feedback gain g for the completely mixed 
state pn with N = 100. One can clearly see the variances 
change symmetrically about the feedback gain. Figure 
[8][b) also shows the noise reduction due to the feedback 
for pioo- Here we define the noise reduction parameter 
of the unpolarized spins as £ unp = V(mi)/V(mo), which 
corresponds to the unconditional squeezing parameter of 
the polarized spins. It is suggested that the optimal gain 
g' is ±l/vJ in single feedback for pioo, and the result is 
consistent with the above estimation of J ~ 0(1/ y/J). 



